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Abstract: Supersymmetric extensions of the Standard Model predict a large en- 
hancement of the Wilson coefficients of the dimension-five electromagnetic penguin 
operators affecting the direct CP violation in Ki — >■ 7r°e"'"e~ and the charge asymme- 
try in — > 7r^£^£~. Here we compute the relevant matrix elements in the chiral 
quark model and compare these with the ones given by lattice calculations. 
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1. Introduction 

Rare kaon decays provide an ideal place both to test the Standard Model (SM) and to 
unravel new physics beyond it [|I| . The origin of CP violation is still an open question 
in modern particle physics. Dimension-five operators including the electromagnetic 
and chromomagnetic penguin operators (EMO and CMO) play important roles in 
these studies since the CP- violating effects from these operators are suppressed in the 
SM but could be enhanced in its extensions 0, P). In fact present experiments, 

HyperCP and KLOE and planned ones, NA48b pl, are going to substantially 
improve the present limits on the Wilson coefficients of these operators by studying 
CP-violating asymmetries in Svr, -k-k^ and in ir^ii {i = e, /i). 

As we shall see, although it is hard to test the SM now it is possible to probe 
interesting new physics scenarios. To this purpose it is necessary to know hadronic 
matrix elements accurately: we address this issue in a particular bosonization scheme. 

The weak effective Hamiltonian, contributed by EMO and CMO, can be written 
as B§ 

Kff = C7+(/.)Q+(/i) + C-{f,)Q-{fi) + C;(/i)Q+(/i) + C;{f,)Q-{fi) + h.c, (1.1) 
where ^ are the Wilson coefficients and 

Ql = ^{^lo^Ar ± SRa^,dj:)F^\ (1.2) 
Ibvr 

= -^(sLa^i.taC?i? ± SE.O^,ytadj}iG^^ . (1.3) 
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Here a^i, = z/2[7^,7j,]. The SM structure, SU{2)l x f/(l), imposes a chiral suppres- 



sion for the following operators |10|, |TT 



njSM T/* 
'^efT = —?=VtdVts 



V2 



9 



+Ci2^{mdSL(y^,ydR + msSRa^,ydL)F^'' +h.c., (1.4) 



and 



Cuim^) = ^TY^^ln X - , , (1.5) 



^ , , x2(2-3x), 8x^ + 5x''-7x ,^ ^. 

Ci2(mtyj = — ; :i-ln a; ; — , (1.6) 



where x = rn^/m^ and ta are the 5'f/(3)-matrices. However, as we shall see, new 
flavour structures in the supersymmetry-breaking terms allow us to avoid the chiral 



suppression for the operators in eq. (|1.4|) 



Among rare kaon decays, the flavour- changing neutral current (FCNC) transi- 
tions K — > 7r£"'"£~, induced at the one-loop level in the SM, are well suited to explore 
its quantum structure and extensions |]12|, |13|, IK. The decay — >• 7r°e+e~ receives 
contributions from three sources [jl], direct CP violation, indirect CP viola- 

tion due to K^-K^ mixing, and CP conservation from the two-photon rescattering 
in Kl — > '7r''77. Therefore, once long-distance effects have been carefully disentan- 
gled new physics, induced by the operators in eq. ( p..2| ), can be probed in this 



channel. Analogously the charge asymmetry in n^i^i^ could be enhanced 

by a large Wilson coefficient of the operator in eq. ( |1.2| ) [|T6|. Recently it has been 
shown that also T-odd correlations in charged K;4-decays depend upon the effective 
Hamiltonian in (|TTTp |T^. 

We thus consider here the matrix element {n'^\Q^\K'^) to determine the observ- 
ables discussed above. In order to evaluate the bosonization of the EMO we exploit 
the chiral quark model, which provides an effective link between QCD and low en- 
ergy chiral perturbation theory. This is particularly interesting since the first lattice 
calculation of the matrix element {7r^\Q^\K^) has been done in Ref. and thus a 
comparison of the two methods can be performed. This might be useful in general 
to understand the extent of validity of the two approaches in the evaluation of other 
matrix elements such as the penguin operator. 



2. The chiral quark model 



The chiral quark model (xQM ) [0 has been extensively used to study low energy 
hadronic physics involving strong and weak interactions [|19|, |21|, |22|, ^ ^ 



Note that the interactions among mesons proceeds in this model only by means of 
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quark loops: starting from the short- distance effective Hamiltonian in terms of quark 
operators (such as four-quark operators, EMO, and CMO), the xQM allows us to 
deduce the low energy effective lagrangian in terms of the input parameters of the 
model. 

In the xQM H^, a term that represents the coupling between the light (con- 
stituent) quarks and the Goldstone mesons 



MQiqnUqL + qLU'-qR) 



(2.1) 



has been introduced into the QCD lagrangian. The Goldstone meson fields, 4>{x), 
are collected in a unitary 3x3 matrix U = exp(i//^A ■ 0(x)) (where the A'^'s are the 
3x3 Gell-Mann matrices and ~ 93 MeV) with detf/=l, which transforms as 



U ^ VnUV^ 

under chiral SU(3)lxSU(3)r transformations (Vl, Vr), and 



(2.2) 



V2 



A • 0(x) 



TT 



TT 



2lL -L HI. 



[2.3) 



In the presence of the term (|2.1| ), it is convenient to use new quark fields, Ql 
and Qr, called "rotated basis" , defined as follows 



Ql = iqi, 
Qr = ^^Qr, 



Ql = qit 
Qr = Qr^, 



with ^ chosen such that 



u=e- 

The chiral SU(3)i^xSU(3)/j transformation 

e(x) ^ VRCix)h+{x) = h{x)C{x)V+ 



(2.4) 
(2.5) 

(2.6) 



defines the compensating SU(3)y transformation h{(j){x)), which is the wanted in- 
gredient for a non-linear representation of the chiral group. Then Ql, ij's transform 
as 

Ql ^ h{x)QL, Qr ^ h{x)QR, (2.7) 

while the term (p. 11) 



-MgiqRUqL + qiU^qR) = -Mq{QrQl + QlQr) 



(2.8) 



3 



is invariant. Therefore, the quark fields Ql, r can be interpreted as "constituent 
chiral quarks" and Mq as a "constituent quark mass" . 

Now in order to evaluate the bosonization of the EMO, we firstly write down the 
EMO using the "rotated basis" in the Euclidean space 

+Q (^^^A+em. + ^^e+A+e+m,^ a,,Q C*^moF"', (2.9) 
where Xij = 6i36j2, and 

G 6 

Cemo = ^^TT^At C*i2, A( = VtdV^l. (2.10) 
V2 Svr^ 

Here we use the form of EMO in the SM [eq. ( |1.4| )]. It is very easy to extend it to 
the general form in eq. ( |1.1[ ). 

Then the effective action induced by the EMO can be written as follows 

r^(AM) = -i I d'xTr j^'^j ^exp [-r(p| + Mj)] exp(-rD'), (2.11) 

where Tr is the trace over colour, fiavour and Lorentz space, V is defined in ( |A.25| ), 



and the integral over r is introduced by using the proper time method The 



detailed derivation for eq. ( p.ll|) has been shown in the Appendix, and dimensional 



regularization has been used for the involved divergences. Expanding exp(— rP') 
in powers of r, and integrating over the momenta, one can get the effective action 
in powers of r, and the corresponding coefficients are the so-called Seeley-DeWitt 
coefficients. Then the effective lagrangian can be obtained by integrating out r. The 



standard procedure can be found in Refs. |2^, 0. If we set Fi = F2 = J,_iu = in 
V [see ( |A.25| ) in Appendix], which implies that the EMO is switched off, eq. ( ^.llj ) 
will give the same effective lagrangian as in Ref. |T^. Here we are concerned about 



the effective lagrangian generated from the EMO, which is relevant to — 7ri~^i 
transitions. Thus at the leading order we get 



-^^Imo = '^^^"^ CEMoiniAUL^L, + ms\L^L,U^)F>''' + h.c, (2.12) 



where = iU^D^JJ ^ Nq is the number of colours, and {A) denotes the trace of A in 
the flavour space. Likewise, the corresponding effective lagrangian from the general 
form of the EMO in eq. (|1 . 1| ) is 



^C^MO = ^V^TfiC7±(Af/L^L, ± \L,UU^)F^^ + h.c. (2.13) 
Svr lovr 

where C^y^Q (-^emo) generates parity-even (odd) transitions. 
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The matrix elements of the EMO between a and a 7r° can be written as 

{AQ^\K') = ^^^|^«F..5t, (2.14) 
' Ib-K niK 

{AQ-,\K') = Q. (2.15) 
Then from eqs. ( [1.1| ) and ( |2.13| ), we can obtain 

Setting Mq =0.3 GeV, we have Bt = 1.31, which is consistent with Bt = 1.18±0.09 
found in the lattice and Bt ~ 1 in Ref. [|1^], and the range \Bt\ = 0.5 ~ 2 



adopted in Ref. [^. Our theoretical error on B^ in (|2.16| ) has two sources: i) from 
the quark mass Mq, which we believe it is very small, ~ 10%, and ii) from higher 
order corrections in the xQM, generated by large- iVc gluonic interactions. We have 



evaluated this contribution using the standard techniques in Refs. |]T9|, ^ finding 
the correction to ( f^.l6D 



9Nr, MA 



(2.17) 



The size of the gluon condensate cannot be simply related to the one which appears 
in the QCD sum rule |2^. However terms like the one in ( CT ), but with lar ger 



coefficients, correct also the leading order predictions for the Lj's and [|T9|. Model 
consistency and the phenomenologically successful predictions of the leading order 
evaluation, lead us to the reasonable expectation that the gluon correction in ( p.lTj ) 



cannot exceed ~ 30% and so consequently we can very conservatively estimate the 
error in this way on B^, i.e. Bt = 1.31 ± 0.4. 

We stress that the agreement with the lattice is found for natural values of the 
chiral quark model. So we can be quite confident in this result. 

3. K ^ 7r£+£- 

The decay width of Kl 7c^e~^e~ induced by the EMO is given by 

Bt{Kl vr°e+e-)EMO = 8.9 x 10^ GeV^ B^ llmC+p. (3.1) 

To obtain an interesting bound on ImC^ we improve our error on Bt by considering 
also the lattice results 0. Thus from the experimental upper bound |^ 

Bt{Kl 7r°e+e") < 5.1 x 10~^°, (3.2) 

we get 

|ImC+| < 1.8 X 10-^ GeV"^ (3.3) 
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at 80% C.L. 



It is known that n^i^i is dominated by long-distance, charge-symmetric, 

one-photon exchange [[12| , [30|1 . This piece can be written as 



TT 



m^(47r) 



rW+(2;)(p^ +p,)^M(p_)7^t;(p+), (3.4) 



where z = {pk — P-kY /^'k-, and the general form factor Wj^{z) has been shown in 
Ref. [^. The piece induced by the EMO will interfere with the imaginary part of 
Wj^{z), which arises from the two-pion intermediate state The asymmetry is 

then written as 



2T I 



EMO 



\V{K+ 



TT 



T{K- 



TX 



-tn\ 



EMO 



V{K- 



TT 



V{K- 



TT 



(3.5) 



Interestingly, with a kinematical cut z > Am^/mj^, the charge asymmetry in eq. 
( p.5|) could be substantially enhanced ||16||. Thus from eqs. ( |3.4| ) and ( p.5|) , and using 
the upper bound of |ImC^| given in eq. ( p.3| ), we can find the charge asymmetry for 



e, /i as 



r-p N EMO 



EMO 



(5r\ 

— < 4.5 X 10" 
2T J 

/ u 



(3.6) 



without the kinematical cut for z, and 



2f < 1-2x10- 



— < 1.3 X 10--* 
2T J 



(3.7) 



with the cut z > Aml/mj^. Note that, differently from Ref. |T6|, here we only use the 
experimental bound of Bt{Kl 7T^e~^e~) to estimate the charge asymmetry in both 
electron and muon mode. So we are neglecting possible lepton-family violations. 



4. Limits on new flavour structures 

From eq. (|2.12| ), one can get ImC^ in the SM 



ImC+|S^ = —^{m, + m,)\lmXt Cu\. (4.1) 



V2 



Due to the smallness of ImXt ~ 10"'^, this contribution from the SM is strongly 
suppressed, and far smaller than the upper bound (|3.3|) . Therefore in the following 
we turn our attention to physics beyond the SM. 

Among the possible new physics scenarios, low energy supersymmetry (SUSY) 
represents one of the most interesting and consistent extensions of the SM. In 



generic supersymmetric models, the large number of new particles carrying flavour 
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quantum numbers would naturally lead to large effects in CP violation and FCNC 
amplitudes [^]. Particularly, one can generate the enhancement of g at one- loop, 
via intermediate squarks and gluinos, which is due both to the strong coupling con- 
stant and to the removal of chirality suppression present in the SM. Full expressions 
for the Wilson coefficients generated by gluino exchange at the SUSY scale can be 
found in Ref. [Q. We are interested here only in the contributions proportional to 
rrig, which are given by 

C^svsYi^s) = ^^^^ mR)2i ± {SEkTu] Fsvsy{x,,), (4.2) 
rrig 

= ^^^^ [(^lk)2i ± (5?r)i2] Gsusy(x,,), (4.3) 



rrig 



where (^lr)^^ = {M^)i^j^/m'^ denotes the off-diagonal entries of the (down-type) 
squark mass matrix in the super-CKM basis, Xgg = rn^/m?^ with rrig being the average 
gluino mass and rrtq the average squark mass. The explicit expressions of -FsusY(a;) 
and GsusY(a;) are given in Ref. 0, but noting that they do not depend strongly on 
X, it is sufficient, for our purposes, to approximate -Fsusy(3^) ~ -Psusy(I) = 2/9 and 
Gsvsy{,x) ~ Gsusy(I) = —5/18. In any case it will be easy to extend the numerology 
once Xgq is better known. Also the determination of the Wilson coefficients in eqs. 
( [4.2| ) and ( ^.3| ) can be improved by the renormalization group analysis [0, ^. Then 
by taking rrig = 500 GeV, = 174 GeV, rrib = 5 GeV, and fj, = rric = 1.25 GeV, we 
will have 

|lmC+|'^''' = 2.4 X 10-^GeV-i \M{S^k)2i + (0*i2]| • (4-4) 
From eq. (|3.3| ), we obtain 

|Im[(5°I,)2l + (5?R)y|<7.7xlO-^ (4.5) 

comparable with the one given by the lattice calculation B. 



5. Conclusions 

To conclude, supersymmetric extensions of the SM may enhance the Wilson coef- 
ficients of the electromagnetic penguin operators. This leads to interesting phe- 
nomenology to be studied: the direct CP violation in Kl n^e'^e' and the charge 
asymmetry in iT^i~^i~ . To this purpose we evaluate the relevant matrix ele- 

ment in the xQM. Interestingly we find a very good agreement with lattice results 
for the natural parameters of the model The present experimental upper bound 
of Bt{Kl — > 7r°e+e~) allows to obtain an upper bound of |ImC+|, and thus to predict 
the upper bound of the charge asymmetry in — > 7i^i~^i~ induced by EMO. The 
analysis shows that the predictions for the relevant matrix elements are solid and thus 
high precision measurements of CP-observables might probe interesting extensions 
of the SM. 
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A. Appendix 



Here we present the derivation for eq. in the xQM. Including the constituent 



quark mass term in eq. the strong lagrangian in the rotated basis [eq.(p.4|)] 

and in the Euclidean space is (after we switch off contributions from the EMO) 



where is the gluon fields strength tensor, and De the Euclidean Dirac operator 

De = i^N^ + M = -i^{d^ + A^) + M, (A.2) 

with 

= iG, + M = - 75A) - Mq. (A.3) 



Note that, in the present paper, we use the same notations as in Ref. |2y] and so 
for the Euclidean quantities, 7+ = 7^, {7;,,7i.} = 25^^, and a^y = -i/2[7^, 7^,]. The 
external vector and axial-vector fields now appear in and 

r, = lieid, - I r,)i + i{d, - I (A.4) 

= t[eid, - 1 - m - ^ (A.5) 

and 

E = eMe+cMt A = eMe-cMt (a.6) 

Here A4 is the current quark mass matrix, and 

r;r = -r;., C = eM, s+ = s, a+ = -a, m^ = m. (a.7) 

The E- and A-terms break chiral symmetry explicitly. 

The Euclidean effective action WE{U,r,l, Jli, Mq) is obtained as follows 



expWE{U,r,l,M,MQ) = VG^exp (^-Ig^I.G^^^ exprs(AM), 



(A. 



where Z is the normalization factor, and 



expri5(^,M) = j VQVQexp J (I^xQDeQ = detDE- (A.9) 
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Since we are concerned with the non-anomalous part of the effective action, we have 

Te{A, M) = detD-^DE, (A. 10) 



with 



Using the technique of the heat kernel expansion one can derive the effective 
strong lagrangian starting from ( [A. 101) , which has been discussed extensively in the 



hterature. 

Now we switch on the EMO. Note that this operator has been expressed using 
the rotated basis in eq. (|2.9| ); it is thus easy to know that ( |A.10| ) should become 



with 



and 



Ve{A, M) = detD+'DE', (A. 12) 



De' = De + J, D+' = D+ + J+, (A.13) 
J '^fii/Jfiuj J '-^lii'J^ui (A- 14) 



- (^^^A+^m. + i^e+A+e+"^d) Qmo^.., (A.15) 

J^u = (75 -75). (A. 16) 
Thus, one can get 

D+'De' -Ml = -V^V^ + E + Fi + F2, (A.17) 

with 

E = iMQ-i^-i^i^ - '-a^,R^,, (A. 18) 

Fl = -'Jf.aalBdf.JalB + -It^C^a^ilb^fi, Ja^} " Mq(T^^{J^u + J+^), (A. 19) 

F2 = -4z7/.-fM-V„ (A.20) 

and 

R,. = iG^, - I (l±^i+FR^,i + l^^FL^^^e^ ■ (A.21) 

Here we set S = A = 0, c?^ is the covariant derivative with respect to the F^- 
connection, i.e. d^j_A = d^A + [F^, A], and the relation 

[7m, o-q/?] = 2i(5^a7/3 - 5^/37a) (A.22) 
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has been used. We only include the linear terms of J^u in Fi and F2 because we are 
concerned about the 0{Gf) AS* = 1 transitions. 

Starting from ( |A.12| ), and in terms of the proper time method we have 



1 r riT 
Te{AM) = -- d^xTi / —{x\eM-rDtDE')\x), 



(A.23) 



where the trace is taken in colour, flavour, and Lorentz space. By inserting a complete 
set of plane waves and using ( |A.17| ), we obtain 



Te{AM) = -^J d'xTr l^^J ^^^P [-ripl + M^Q)] exp(-rl?'), (A.24) 
V' = E-V-V + F,- 2ipE ■ V + F2 + A^^pEJ^.u. (A.25) 



where 
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